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Abstract. We present a statistical study of the presence 
of galaxy warps in the Hubble deep fields. Among a com- 
plete sample of 45 edge-on galaxies above a diameter of 
l."3, we find 5 galaxies to be certainly warped and 6 
galaxies as good candidates. In addition, 4 galaxies reveal 
a characteristic U-warp. Compared to statistical studies 
of local warps, and taking into account the strong bias 
against observing the outer parts of galaxies at high red- 
shift, these numbers point towards a very high frequency 
of warps at z ~ 1: almost all galaxy discs might be warped. 
Furthermore, the amplitude of warps are stronger than for 
local warps. This is easily interpreted in terms of higher 
galaxy interactions and matter accretion in the past. This 
result supports these two mechanisms as the best candi- 
dates for the origin of early warps. The mean observed 
axis ratio of our sample of edge-on galaxies is significantly 
larger in the high-z sample than is found for samples of 
local spiral galaxies. This might be due to disk thickening 
due to more frequent galaxy interactions. 

1. Introduction 

Warping of outer galactic disks is a very common phe- 
nomenon, with amplitude increasing with radius. Warps 
are mostly conspicuous in the HI plane, since the atomic 
gas extends up to lOOkpc, and its radius is often several 
times that of the optical disk (e.g. Bosma 1978, Briggs 
1990). However warps start at the outskirts of optical 
disks, and careful investigation of edge-on galaxies have 
revealed that a significant percentage (40-50%) of stel- 
lar disks are also warped (Sanchez-Saavedra et al. 1990, 
Reshetnikov & Combes 1998). 

The origin of warps has been a puzzle for a long time, 
since orbits slightly inclined with respect to the plane pro- 
cess around the disk rotation axis, and the precession is 
differential: the life-time of warps would be much less than 
the Hubble time. Several mechanisms to trigger and main- 
tain warps have been proposed, including tidal interac- 
tions between galaxies, inter-galactic magnetic fields, disk- 
Send offprint requests to: V.P. Reshetnikov resh@astro.spbu.ru 



halo interaction and matter accretion (see the review by 
Binney, 1992). The forcing of a misaligned halo does not 
appear to be a promising solution, since the halo and disks 
align over a short time-scale (Nelson & Tremaine 1995, 
Binney et al. 1998). On the contrary, galaxy interactions 
have been revived as an explanation, since even the proto- 
typical isolated warped galaxy NGC 5907 has been found 
to interact with small companions (Shang et al. 1998). 
The other likely alternative is gas accretion, since it has 
been realized that galaxies are still significantly growing 
in mass during all their lives, and are surrounded by a net- 
work of filaments accounting for the Lyman-a absorbers 
(e.g. Dave et al. 1999). The outer parts of galaxies could 
be constantly accreting gas at an inclined angular momen- 
tum, explaining the persistence of warps. The dark matter 
is also constantly raining down onto the dark galactic halo, 
with a skewed angular momemtum, making the galaxy re- 
orient completely in a few Gyrs (Jiang & Binney 1999). 

If these two mechanisms, matter accretion and tidal in- 
teraction with companions, are the main drivers of warps, 
it is very likely that warps were even stronger and more 
frequent in the past. Knowing this frequency at high red- 
shift would help to discriminate between models of the 
warp origin. We thus undertake such an investigation in 
the Hubble Deep Fields. The main goal of our work is to 
study the frequency of optical warps at z ~ 1. At present, 
HI-21cm observations are not sensitive enough (and have 
insufficient spatial resolution) to probe high redshift galax- 
ies. 

Throughout this paper, we adopt a flat cosmology with 
qo = 0.5 and Hq = 70 kms _1 Mpc _1 . If not differently 
specified, magnitudes are expressed in the AB system (Oke 
1974). 

2. The sample 

The Hubble deep fields north and south (HDF-N and 
HDF-S, respectively) are the best available sources of data 
on the distant Universe (see Ferguson et al. 2000 for a com- 
prehensive review). The two fields give the possibility to 
study distant field galaxies to the faintest levels currently 
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achievable. Therefore, we selected HDF-N and HDF-S to 
study optical warps in remote galaxies. 

2.1. Selection 

We select our sample of distant edge-on galaxies on the 
basis of visual inspection of original HDF-N and HDF-S 
frames in the F814W filter (hereafter Isia)- Visual selec- 
tion is probably the best way to recognize true strongly 
inclined spirals among a huge number of asymmetric and 
peculiar objects in the two crowded deep fields. 

The resulting sample of selected objects consisted of 
50 edge-on galaxies. Tables 1 and 2 summarize the main 
characteristics of the galaxies. 

The columns of Table 1 arc: number according to the 
catalogue of Fernandez-Soto et al. (1999); coordinates (in 
pixels); Williams et al. (1996) identification; 7gi4 magni- 
tude in the AB photometric system; spectroscopic redshift 
(Cohen et al. 2000; Hogg et al. 2000); photometric redshift 
(Fernandez-Soto et al. 1999); photometric redshift accord- 
ing to Fontana et al. (2000); best-fit spectral type of galaxy 
(Fernandez-Soto et al. 1999); index of warp (see sect. 3. 2). 

The columns of Table 2 are: the name according 
to the Stony Brook catalogue (Fernandez-Soto al. 1999) 
(the numbers mean coordinates in pixels); I%u magni- 
tude in the AB photometric system; photometric redshift 
(Fernandez-Soto et al. 1999); photometric redshift accord- 
ing to Fontana et al. (2000); best-fit spectral type of galaxy 
(Fernandez-Soto et al., 1999); index of warp. 

2.2. Completeness 

Visual classification is possible for relatively large galaxies 
only. Thus one can expect that our sample is a diameter- 
limited sample. In order to estimate the completeness limit 
we applied the well-known V/V ma x test (Schmidt 1968). 
For a diameter-limited sample the ratio between the vol- 
ume of the sphere determined by the object distance (V) 
and the maximum volume of the sphere in which the ob- 
ject could lie and still have been included in the sample 
(Vmax) is V/V marE = (do/d) 3 , where d is the apparent an- 
gular diameter of the galaxy and do is the minimum diam- 
eter of the sample selection criterion (e.g. Thuan & Seitzer 
1979). For objects distributed uniformly in space the av- 
erage Value Of V/V max should be 0.5+1/V12 x N, where 
N is the number of sample galaxies. 

We determined the apparent major axis angular di- 
ameter (d) within the isophote f^(Isu) = 26.0 for all 
the sample galaxies. Then, changing do and comput- 
ing the V/V max ratio, we found that for d =l."33 
(V/V ma:E )=0. 50+0.04. Therefore, our sample is nearly 
complete for the galaxies with diameter d>l."3. (Note 
that we used the isophotal sizes of galaxies for the V/Vmax 
test. This is not quite correct due to the influence of the 
cosmological dimming and if-correction on the apparent 
sizes of galaxies. Note also that surface brightness evolu- 



tion - e.g. Lilly et al. 1998 - must partially compensate 
the above mentioned effects.) 

Only five galaxies (nll2, n246 in the HDF-N and 
s0932, s2504, s2816 in the HDF-S) show major axis di- 
ameters smaller than l."3. We will consider the rest of 
the sample - 45 galaxies - as a complete sample. 

2.3. General characteristics 

Fig.l shows the distribution of the selected objects within 
the HDF-N and HDF-S. Uniform angular distribution of 
galaxies is evident in the figure. HDF-N contains more 
edge-on galaxies (32) than HDF-S (18). This difference 
probably reflects the general excess of galaxies in the 
northern field at z < 1 (Gwyn 1999). 

We present available spectroscopic and photometric 
redshifts of the galaxies in tables 1 and 2. Unfortunately, 
we have spectroscopic redshifts for only about half of the 
objects (14) in the HDF-N. For the remaining galaxies we 
will use the mean photometric redshifts from the cata- 
logues of Fernandez-Soto et al. (1999) and Fontana et al. 
(2000). The edge-on orientation and, therefore, strong in- 
ternal absorption can distort the observed spectral energy 
distribution of a galaxy and can introduce some system- 
atic effects in the derived photometric redshift. To check 
for this possible effect, we compared two catalogues of 
photometric redshifts with spectroscopic data. The two 
sets of photometric redshifts demonstrate relatively good 
agreement (Fig. 2): 

{zfsoto - Zf on t) = +0.00 ± 0.31(a) (7V=50) or 
(zfsoto - z font ) = +0.05 ± 0.13 (z < 1.8, 7V=46). 
The comparison of spectroscopic and photometric red- 
shifts does not show any systematic effects (Fig. 2): 
(z - z fsoto ) = +0.02 ± 0.26 {N=U), 
(z - z }ont ) = -0.03 ± 0.31 (JV=14). 

The largest differences between photometric and spectro- 
scopic z are for n716 and n727. These two galaxies are 
very close in projection (2. "9 or 17 kpc at z = 0.92) and 
show close spectroscopic z (see table 1). We suppose that 
these galaxies form a real physical pair and that their spec- 
tral energy distributions are strongly distorted by tidally- 
induced star formation. 

Fig. 3 presents the redshift distribution of edge-on 
galaxies. The distribution is peaked at z ~ 0.8 — 0.9 (ta- 
ble 3) and is located in the low-redshift wing of a global 
redshift distribution of galaxies in the deep fields (bottom 
part of Fig. 3). This is probably explained by our diameter- 
limited selection of the sample. 

Fig. 4 shows the distribution of apparent axial ratios 
- b/a - of the sample objects. The b/a values were es- 
timated from the ellipse-fitting of outer contours of the 
galaxies. Our b/a estimates are very close to those pub- 
lished in the Williams et al. (1996) catalogue. The mean 
rat io fe / / a( P r ^; work , ) is 1.01±0.21(<r). The average value of 

b/ a (Williams) x - ' & 

axial ratio - 0.32 (table 3) - significantly exceeds the b/a 
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N 


Mosaic (X,Y) 


Williams et al. 
(1996) ID 


^814 


z 


Zfsoto 


Zfont 


Type 


Warp 


15 


733.9,271.3 


3-923.0 


24.70 




0.40 


0.49 


Scd 


6 


16 


1342.8,274.7 


3-925.0 


24.93 




0.92 


0.88 


Irr 


U 


55 


3893.6,380.4 


4-976.2 


23.32 




0.20 


0.15 


Irr 




81 


3232.6,433.0 


4-618.0 


23.94 




0.64 


0.60 


Irr 




112 


2764.7,508.1 


4-386.0 


26.09 




1.08 


0.87 


Irr 




120 


1368.7,534.7 


3-786.0 


24.27 




1.60 


1.60 


Irr 




153 


673.6,598.0 


3-761.1 


24.94 




0.76 


0.69 


Irr 




246 


2154.6,809.8 


4-70.0 


26.99 




1.48 


1.69 


Irr 




273 


2107.1,872.3 


4-33.0 


22.47 


0.905 


0.68 


0.80 


Ell 




304 


1747.1,951.1 


3-593.1 


25.60 




1.760 


1.79 


Scd 


u 


319 


2178.7,989.2 


4-85.1 


24.41 


0.961 


0.92 


0.90 


Irr 


u 


450 


929.2,1329.3 


3-398.1 


24.25 




0.72 


0.67 


Irr 




476 


2430.9,1390.5 


4-232.12 


22.77 


0.421 


0.48 


0.45 


Irr 




506 


1899.5,1472.4 


3-331.0 


24.47 


0.751 


1.08 


0.75 


Irr 




534 


3812.9,1545.4 


4-926.1 


23.82 




1.04 


0.85 


Sbc 


5 


632 


3483.8,1771.9 


4-745.0 


25.80 




1.60 


1.88 


Irr 




671 


2166.4,1877.4 


4-89.0 


23.93 


0.681 


0.76 


0.65 


Irr 




716 


3849.0,2034.0 


4-946.0 


22.61 


0.944 


0.68 


0.45 


Sbc 


5 


727 


3792.2,2078.6 


4-916.0 


23.93 


0.904 


0.16 


1.88 


Scd 


6 


733 


3628.1,2104.3 


4-829.1 


25.04 




0.96 


0.91 


Irr 




749 


1458.6,2208.5 


2-339.1 


22.60 


0.851 


0.88 


0.95 


Sbc 


U 


774 


2521.0,2361.5 


1-34.1 


21.22 


0.485 


0.68 


0.40 


Sbc 




805 


172.2,2492.4 


2-1022.0 


25.33 




0.72 


0.59 


Irr 




817 


2793.4,2550.5 


1-57.11111 


22.19 


0.485 


0.68 


0.50 


Sbc 


4 


886 


1835.2,2838.4 


2-133.0 


25.25 




1.00 


0.94 


Irr 




888 


196.0,2852.5 


2-1018.0 


23.70 


0.559 


0.64 


0.56 


Irr 




898 


1958.2,2909.9 


2-65.1 


25.89 




0.92 


0.95 


Irr 


6 


899 


181.2,2923.1 


2-1023.1 


22.94 


0.564 


0.68 


0.60 


Scd 




938 


765.0,3128.9 


2-702.1 


23.03 


0.557 


0.52 


0.55 


Scd 




979 


327.1,3397.0 


2-950.0 


23.45 


0.517 


0.44 


0.50 


Scd 


4 


1027 


1230.6,3650.7 


2-432.0 


23.88 




0.92 


0.87 


Irr 




1031 


1590.7,3680.9 


2-270.2 


25.08 




1.20 


0.93 


Scd 





Table 2. Edge-on galaxies in the HDF-S 



Name 


^814 


Zfsoto 


Zfont 


Type 


Warp 


SB-WF-0501-0760 


25.07 


1.11 


1.31 


Irr 




SB-WF-0564-0515 


25.27 


1.01 


1.03 


Irr 




SB-WF-0566-1028 


25.67 


1.79 


2.52 


Scd 




SB-WF-0578-2216 


23.16 


1.08 


0.88 


Scd 


6 


SB-WF-0747-1929 


25.31 


1.07 


0.88 


Irr 




SB-WF-0932-2115 


25.58 


1.03 


0.88 


Scd 




SB- WF- 1085-0562 


25.51 


1.07 


1.19 


Irr 




SB- WF- 1404-2925 


23.09 


0.54 


0.45 


Sbc 




SB-WF-2340-1898 


24.21 


0.55 


0.56 


Irr 




SB-WF-2353-1914 


24.36 


0.30 


0.31 


Scd 




SB-WF-2504-3117 


24.80 


1.20 


1.09 


Irr 


4 


SB-WF-2661-4090 


24.10 


0.92 


0.92 


Irr 




SB-WF-2691-2966 


25.05 


0.49 


0.42 


Irr 




SB-WF-2816-2895 


26.04 


0.47 


0.63 


SB1 




SB-WF-3053-1963 


26.07 


1.18 


0.91 


Scd 


5 


SB-WF-3329-2173 


24.54 


0.27 


0.32 


Scd 




SB-WF-3458-1026 


23.40 


0.69 


0.69 


Scd 


6 


SB-WF-3685-2040 


26.63 


2.43 


2.03 


SB1 
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Fig. 1. Distribution of edge-on galaxies within the two fields. Crossed circles - galaxies with possible warps (see tables 
1 and 2). 
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Fig. 2. Left: comparison of photometric redshifts of our sample edge-on galaxies according to Fernandez-Soto et al. 
(1999) and Fontana et al. (2000) catalogues. Right: comparison of spectroscopic and photometric redshifts for 14 
edge-on galaxies. 



ratio for local spiral galaxies. According to Guthrie (1992) 
the true axial ratio of spiral disks is only 0.11-0.13. As- 
suming a very conservative true axial ratio of our sample 
of edge-on galaxies b/a = 0.2, we can estimate that the 
observational value (0.32) corresponds to a rather large 
inclination of the galaxies (i <~ 75°). Such a large depar- 
ture from edge-on orientation must lead to very strong 
asymmetry of the surface brightness distributions along 
the minor axis. However we do not see any such signif- 
icant asymmetries among the sample galaxies. Another 
possible reason is that the sample objects possess a large 
intrinsic axial ratio (b/a w 0.3) - about twice that of local 
spirals. This conclusion is consistent with the facts that 
1) interaction and merger rates were significantly higher 
in the past (see Le Fevre et al. 2000, Reshetnikov 2000 



and references therein) and 2) galaxy interactions lead 
to thickening of the stellar disks (Reshetnikov & Combes 
1997; Schwarzkopf & Dettmar 2000; see also Karachent- 
sev & Karachentseva 1974). Of course, our conclusion is 
extremely tentative since we must check possible influence 
of the PSF on the observed thicknesses of edge-on galaxies 
in the HDF. 

The mean isophotal diameter of the sample galaxies 
within At(/ 8 i4) = 26.0 is 2."15±0."9 or 12.9 kpc at z = 0.9. 

The rest-frame absolute magnitudes in the B passband 
were estimated following Lilly et al. (1995): 
M(B) = J 814 - 51og(D L /10pc) + 2.51og(l + z) + [B - 
hu.z), where Dl is the luminosity distance and B — Isi4, z 
is the -ftf-correction color. For our purposes we used the 
B — hu.z term for Irr galaxy (Lilly et al. 1995). At 
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Fig. 3. Top: rcdshift distribution of edge-on galaxies (solid 
line shows galaxies with photometric redshifts, dashed line 
- with spectroscopic determinations). Bottom: distribu- 
tions of the photometric redshifts in the HDF-N (solid 
line) and HDF-S (dashed line) according to Fernandez- 
Soto et al. (1999) and Fontana et al. (2000) catalogues. 

Table 3. Mean characteristics of the sample 



Parameter 


HDF-N 


HDF-S 


N+S 


Number 


32 


18 


50 


</814> 


24.15±1.30(ct) 


24.88±1.02 


24.41±1.24 




0.85±0.38 


0.95±0.54 


0.89±0.44 


{b/a) 


0.31±0.07 


0.34±0.08 


0.32±0.08 


(M B ) 


-18.1±1.2 


-17.7±1.4 


-17.9±1.3 



high redshifts, this term is small since Isi4 is redshiftcd 
down toward the rest-frame B filter (they coincide at 
z = 0.83). The resulting mean rest-frame absolute mag- 
nitude is M(B) ss — 18 (table 3). After correction for in- 
ternal absorption the absolute magnitude must be w —19. 
Therefore, our sample objects are sub-L* galaxies. 

The mean parameters of the sample are summarized 
in table 3. 



3. Statistics of warped galaxies 




0.2 0.4 0.6 0.8 1 
b/a 

Fig. 4. Apparent axial ratios: HDF-N (solid line) and 
HDF-S (dashed line). 

3.1. Could we recognize local warps at high redshifts? 

Searching for optical warps in the disks of distant galax- 
ies is a difficult enterprise. In local galaxies optical warps 
become detectable in the outer parts of their disks. For 
instance, from a sample of 45 spirals with measureable 
warps, Reshetnikov (1995) estimated that the average 
value of the surface brightness level at which optical warp- 
ing becomes detectable is n(R) = 23.0±0.8 (in the Cousins 
R passband) or /i(B) ~ 24. (After correction for the line- 
of-sight integration, this level corresponds to fi(B) w 25 
for face-on galaxy.) Moving such a typical object at z = 0.9 
(average redshift of our sample), this level will fall to 
/x(l8i4) ~ 26 (see Lilly et al. 1995 and note that the K- 
correction term is close to zero at this redshift). This sur- 
face brightness corresponds to ~ 1.5(7 limiting isophote of 
the deep fields and, therefore, we cannot recognize such a 
typical optical warp at z <~ 1. 

In order to directly illustrate this conclusion we have 
redshiftcd the B band image of the local Sb spiral 
ESO 235-G53, which has a remarkable optical warp. The 
image of the galaxy was taken by R. de Grijs with the 
Danish 1.54-m telescope at La Silla, Chile, in July 1994. 
The prominent warp starts at the radius w (1.2 — 1.7) ft 
and reaches s=s 10° near the end of the disk (de Grijs 1997). 
Using the method described in Giavalisco et al. (1996), we 
placed ESO 235-G53 at z = 0.83 where redshiftcd B cor- 
responds to Isu- The result of the redshifting is shown 
in Fig. 5. (Note that this is just an illustrative figure and 
we faded the whole image including the foreground stars.) 
The result demonstrates that even such a strong warp 
cannot be detected. 

We can conclude from the above exercise that among 
distant galaxies we can recognize only extremely strong 
and bright optical warps. Any statistics of warps in the 
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Table 4. Characteristics of optical warps 



Fig. 5. Top: contour map of ESO 235-G53 in the B pass- 
band. The faintest contour corresponds to fJ.(B) = 27.0. 
Bottom: contour map of the same galaxy moved at z = 
0.83. The faintest contour corresponds to = 26.5 

or <~ lcr limiting isophote of the HDF. Isophote step is 
0. m 75. Both axes arc in arcseconds. 



deep fields will give a strict lower limit of the true fraction 
due to a huge observational bias against finding them. 

3. 2. Selection of warped galaxies 

The original 7 8 i4-band images and contour maps of all the 
sample galaxies were examined independently by three of 
us (E.B., F.C. and V.R.). On the basis of eyeball classifi- 
cation, they assigned an " index of warp" (IW) for each ob- 
ject: - no warp, 2 - a certain integral-shaped warp, and 
1 - intermediate case. We also marked U-shaped warps, 
when the plane bends towards the same direction on both 
sides (Reshetnikov & Combes 1998) like the letter "U". 
Then, we summed the personal indices for each object 
and decided to consider those galaxies with total IW=6 
(index 2 for each observer) as true warps and IW=4-5 as 
good candidates for warped disks. There are 5 galaxies 
with true warps, and 6 good candidates in our sample of 
50. Our classification for the galaxies with IW>4 is pre- 
sented in the last columns of tables 1 and 2. Contour maps 
of all galaxies with certain or possible warps are displayed 
in Fig. 6 and 7. Note that these are just illustrative contour 
maps and they do not give the best impression of possible 
warps. A detailed inspection must be made on the original 
pictures. 

These 11 galaxies were processed with a special soft- 
ware, written to automatically detect and measure some 
quantitative characteristics of warps. The software works 
essentially as the WIG program described by Jimenez- 
Vicente et al. (1997). First the images are rotated so the 
galaxies are roughly in a horizontal position and the pixels 
with a S/N>3 are selected. Then the warp curve (mean 



Galaxy 


i> 


Rwarp 


R^warp / h 




(°) 


(arcsec) 




nl5 


4.5 


0.58 


1.5 


n716 
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deviation from the equatorial plane for each radius) is cal- 
culated by fitting a gaussian to each column. Outliers are 
removed using S/N and continuity criteria. The flat part 
of the galaxy is defined by the average position angle of 
the inner FWHM of the galaxy disk. 

The software detects a warp if there are more than 
five consecutive pixels with deviations larger than twice 
the standard deviation in the flat part. In this case the 
first of those pixels is taken as the starting radius of the 
warp. Then, the warp angle is determined for both sides 
(defined as the angle between the equatorial plane and 
the line joining the center with the point with the largest 
deviation from the equatorial plane in the warp curve) and 
averaged to give the average warp angle. 

As a result, we confirmed that 9 of 11 galaxies demon- 
strate integral-like shapes according to our criteria. 

Table 4 summarizes the general characteristics of op- 
tical warps in the sample according to the software. The 
columns are as follows: V ~ mean warp angle, R warp - 
starting point of the warp, and R war p/fr ~~ starting point 
of the warp to exponential scalelength ratio (uncertain 
values are marked by colons). 

3.3. Statistics of warps 

The measured mean warp angle is (ip)—5.°9zk2.°8(o-) or, 
excluding two extreme cases (n716 and s0578), 5.°8±1.°6. 
This value is larger than reported by Reshetnikov (1995) 
or Reshetnikov & Combes (1998) - 3°-4° - for local spi- 
rals. The mean warp angle for distant objects is also 
larger than for local galaxies with strongest optical warps 
- 4.°8±1.°3 (Reshetnikov & Combes 1999). 

The observed fraction of optical warps in the HDF-N 
and HDF-S is 9/50 = 0.18±0.06 (Poisson error). Consid- 
ering the "complete" sample of 45 galaxies, this fraction 
is 9/45 = 0.20±0.07. The observed fraction is certainly an 
underestimate of the true fraction (see sect. 3.1). 

There is a possibility of confusing a warped disk with 
a slightly inclined (from edge-on) disk with a contrasted 
2-spiral pattern (see the simulations in Reshetnikov & 
Combes 1998). This type of confusion should be less im- 
portant at high redshift, since grand-design spirals appear 
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to be rarer for z>0.3 or even absent in the past (cf van 
den Bergh et al. 2000). Possible projection of small back- 
ground objects, for instance for n898 and s0578, can be 
rejected on the basis of smoothness and symmetry of the 
warp curves. However it is possible that in both cases we 
observe accretion of small companions by a large spiral 
galaxy. 

To obtain a rough estimate of the true fraction of dis- 
tant warps, let us assume that warped disks at z <~ 1 obey 
the same distribution of warp angles as the z = spirals - 
ip~ 5 (Reshetnikov & Combes 1998, 1999). Excluding two 
extreme cases with ip = l.°7 and 11. °4, we have 7 galaxies 
with ip = 4.°5 — 9.°2. Assuming that we selected all spirals 
in the two deep fields with ip — 4.°5 — 9.°2, one can infer 
that the true number of galaxies with ip > 2.°5 must be 31 
or -2/3 of all galaxies (2.° 0-2.° 5 is the detection limit of 
optical warps in our survey of local spirals - Reshetnikov 
& Combes 1998). Therefore, the actual fraction of warped 
disks at z — 1 is the same or even larger than at z = 0. 
Given the geometrical factor (warps along the line of sight 
are not detected), practically all distant disks are warped. 

To draw a clearer conclusion about possible warp evo- 
lution with redshift, let us compare the frequencies of 
strongly warped disks at z — and in the deep fields. 
In the two fields we found 8 galaxies with ip > 4.°5 (table 
4). What is the equivalent fraction in the local universe? 
About 40% of spiral galaxies at z = reveal S-shaped 
optical warps with ip > 2° (Reshetnikov & Combes 1998). 
Adopting the ip~ 5 law, one can estimate that the local fre- 
quency of spirals with ip > 4.°5 is 0.7%. Then, using the 
method described in Reshetnikov (2000), we found that 
the volume density of strong warps changes oc (1 + z) m 
with m = 3.8. Assuming Poisson errors for galaxy num- 
bers, we have m = 3. 8^1^ (owing to the uncertainty in 
the local frequency of strong warps, the actual error on 
m may be larger.) Comparing frequencies of galaxies with 
ip > 5.°0, we obtain to = 4.0l°;g. Our data therefore 
demonstrate a fast evolution of the volume frequency of 
strong optical warps at z — 1. The galaxies are therefore 
more frequently warped, and distant warps are brighter 
and stronger. The rate of the evolution is consistent with 
the increase of galaxy interactions and mergers (Le Fevre 
et al. 2000, Reshetnikov 2000), which tends to suggest 
that the source of warping is precisely tidal interaction. 
Schwarzkopf & Dettmar (2001) estimate that tidal per- 
turbations generate the strongest warps. 

As one can find from Tabic 1, the mean pro- 
jected starting point of the stellar warps we detected is 
(Rwar P )=3.7±1.3 kpc or, expressed in exponential scale- 
lengths, 1.4±0.3ft. According to de Grijs (1997), the same 
value for local spirals is 2.1±1.0/i. 

No coherence in the orientation of warps was found. 
This is not surprising, taking into account the large dis- 
tances involved, of the order of gigaparsecs. Galaxies with 
close redshifts were too few to be of any consequence. This 



phenomenon cannot be excluded, however (Battaner et al. 
1991, Zurita & Battancr 1997). 

4. Discussion and conclusions 

The result of the present statistical study is that warps 
were very frequent in the past, and they were of larger 
amplitude. One obvious interpretation comes from galaxy 
interactions that were much more frequent, since they rep- 
resent an efficient trigger of warps. Indeed the frequency 
of galaxy pairs and the frequency of mergers is observed 
to increase as a high power of time, or in (l+z) m , with 
to = 3 — 4 (Le Fevre et al. 2000). This is also observed 
in the fraction of perturbed and distorted galaxies, for in- 
stance in the frequency of tidal tails (Reshetnikov, 2000). 
However, a much lower rate of mergers or mass accre- 
tion was observed for bright and massive galaxies by Carl- 
berg et al. (2000). This might be related to the fact that 
the spheroids or bulge-dominated galaxies appear later in 
redshift, at the expense of disk-dominated and irregular 
galaxies (Kajisawa & Yamada, 2001). This is expected in 
a hierarchical cosmological scenario, where galaxies are 
formed by successive mergers. Galaxies are then more fre- 
quent in the past, per comoving volume. Moreover, due 
to expansion, the density of galaxies was larger (increas- 
ing as (1 + z) 3 for a homogeneous universe). On the other 
hand, galaxies were smaller, and the warps appear to be- 
gin at a smaller radius. At this distance from the center, 
the differential precessing rates of orbits is shorter, and the 
life-time of warps should be smaller (roughly proportional 
to Rwarp, for flat rotation curves). 

Another possible triggering of warps is the action of an 
intergalactic magnetic fields. In a highly conductive ho- 
mogeneous universe, magnetic fields increase as (z + l) 2 
(Battaner & Florido 2001), which means a factor of 4 at 
z = 1. However, at z = 1, with a highly inhomogeneous 
universe the condition of frozen-in magnetic field lines 
should better describe the evolution of magnetic fields. We 
adopt _B/?~ 2 / 3 =constant. Then, for a present-day cluster 
overdensity of 10, the magnetic field strength should be 
2 2 / 3 = 1.6 larger than today. Magnetic fields larger by a 
factor 1.6-4 also constitute a potential explanation of the 
increased warping amplitude in the past. 

The infall of dark and visible matter either on halos 
(Ostrikcr & Binncy 1989) or on disks, producing reori- 
entation of the galactic angular momentum, should have 
been greater in the past. Therefore the infall mechanism 
could also be a dominant effect to explain the high rate of 
z — 1 warps. 

Finally, we would like to stress that our results are ob- 
tained for very small regions of the celestial sphere. More- 
over, the statistics are based on small numbers. Therefore, 
further extended studies of distant edge-on galaxies are 
highly needed to develop our preliminary results. 
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ing Paris, thanks to a one month grant from Paris Observatory. 
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Fig. 6. Contour plots of galaxies with certain or possible warps (see tables 1 and 2). Both axes are in arcseconds (size 
of 1 pixel of the original fields is 0."04). The faintest contour corresponds to n(Isi4) = 26.0, isophotes step is 0. m 44 
(0. m 75 for n817). Uncertain values of z obtained as the mean of two photometric redshifts (see tables 1 and 2) are 
marked as (:). 
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Fig. 7. Contour plots of galaxies with certain or possible warps (continued) . 
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